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Abstract. We have studied the dynamics of quantum correlations such as 
entanglement, Bell-nonlocality and quantum discord between identical as well as 
unidentical atoms interacting with a single-mode cavity field and subject to cavity 
decay. The effect of single atom detuning, cavity decay rate and initial preparation 
of the atoms on the corresponding correlation measures have been investigated. It 
is found that even under strong dissipation, time evolution can create high quantum 
discord while entanglement and Bell nonlocality stay zero for an initially separable 
state. Quantum discord increases while entanglement decreases in a certain time period 
under dissipation for the initial state that both atoms are in the excited state if the 
qubits are identical. For some type of initial states, cavity decay is shown to drive the 
system to a stationary state with high entanglement and quantum discord. 



PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.-a 



Quantum Correlations 



2 



1. Introduction 

Entanglement, nonlocality and quantum discord (QD) are all different but somehow 
related aspects of the quantumness of the correlated composite quantum systems. 
Entanglement is a kind of quantum correlation and determines whether the given 
state is separable or not [JJ. On the other hand, the quantum systems may contain 
quantum correlations other than entanglement. Quantum discord, which is defined as 
the difference between the quantum versions of two analogous classically equivalent 
definitions of mutual information, captures all types of nonclassical correlations 
including entanglement p]. Although entanglement and QD are equivalent for pure 
states [3], QD is found to be nonzero for some separable states [31 H]. Such states were 
shown to be useful in the deterministic quantum computations with one pure qubit 
context [3, E]. Therefore QD is believed a new resource for quantum computation. 
The nonlocality as measured by the violation of Bell's inequalities is also a signature of 
the inseparability of a quantum state and signifies entangled states whose correlations 
cannot be reproduced by a classical local hidden- variable model [7]. However it must 
be emphasized here that the presence of entanglement does not always imply violation 
of Bell's inequalities, while the violation of Bell inequalities implies entanglement. Such 
entangled states violating Bell's inequalities play a central role in some applications in 
quantum information science [8], such as to guarantee the safety of device-independent 
key distrubition protocols in quantum cryptography [9| [TO]. 

The relations among different measures of quantum correlations and their 
dynamical behaviour have been an active area of research with an aim to understand the 
fundamental question of what quantumness is as well as to characterize these correlations 
as useful resources for various quantum tasks P3llI2llI3llIllI3IIElinilIEllISll2nil2IlE2l 
[231 1231 [251 ESI EZl EHl [291 EQl EH [32J [331 [351 [SS] . One of the most important obstacle to 
realizing quantum operations is the interaction between the qubits and the environment 
which tends to wipe out the quantumness of the system pp. Characterization of the 
effect of environment on the system such as decoherence and dissipation is important. 
The dynamics of quantum correlations such as entanglement, quantum discord and Bell 
nonlocality for a two-qubit system subject to various environments have been studied 
by many groups recently E21 H31 OH HSl HSl HZl UEl HS1 1201 I2H 1221 1231 EH I2SI 
I2S1I2Z1I2S1I221I3Q1EI1I321I331E11I351ES]- Entanglement sudden death (ESD) is one 
of the new concepts that arise within the context of these studies [37]. ESD refers 
to death of entanglement between two qubits in a finite time while the coherence of 
the single qubit decays only exponentially. In contrary to entanglement, QD is more 
robust in a decohering environment, actually it was shown that almost all states have 
non-zero quantum discord [38]. QD was shown to decay exponentially in Markovian 
environments, even at finite temperatures [HI [16] and can become zero momentarily 
under non-Markovian dynamics [IB] . QD was also experimentally investigated in an 
all-optical set up recently [39]. Furthermore, QD and entanglement were analyzed for a 
Heisenberg spin chain with quantum phase transitions and QD was found to detect the 
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critical points of the transition at finite temperatures while entanglement could not [f 7J. 
Very recently, Sun et al. investigated the induced quantum correlations between two 
qubits coupled to a common environment modelled by Ising spin chain in a transverse 
field and during the time evolution, non-zero quantum discord is found to be created 
for specific initial states whereas entanglement cannot be induced [12]. In Ref. [21] 
quantum discord is found to be totally unaffected by non-dissipative noise for a long 
time interval under certain conditions, while entanglement decays to zero quickly These 
remarkable properties of quantum discord makes it more practical than the entanglement 
to characterize quantum correlations. 

As mentioned, the nonlocality is curicial for some applications in quantum 
information. Therefore, many theoretical studies have been devoted to the comparision 
of entanglement and the violations of Bell inequalities [211 [3Qj [3TJ [321 EHl [34j [35j [36] . 
It was demonstrated by several authors that the survival time for entanglement is 
much larger than that of the Bell- inequality violation [291 EQl EH EH]. In Ref. |34j, 
the nonlocal entanglement identified by the violations of a Bell inequality is found 
to be preserved during the time evolution of a system consisting of two qubits which 
are embedded independently in zero-temperature bosonic reservoirs and are initially 
in an entangled mixed state. Bellomo et al. showed that Bell inequality might not 
be violated for a state with high entanglement for a two-qubit system subject to 
amplitude damping [35]. Moreover, the nonlocal entanglement between two qubits 
subject to classical dephasing environment modelled as Ornstein-Uhlenbeck noise has 
been investigated in Refs. [301 El] , very recently, and it was demonstrated that the strong 
non-Markovian effect can protect the nonlocal entanglement [301 EI]- Additionally, the 
nonlocal entanglement induced by the qubit-qubit interaction is found to be protected 
against the sudden death despite the strong influence of the environmental noise for 
a certain class of initial states regardless of whether the environment is Markovian or 
non-Markovian [31] . 

Although the dynamical relation between the entanglement and quantum discord 
as well as entanglement and Bell nonlocality for a bipartite system subject to different 
environments have been investigated by a large number of groups as mentioned above, 
further investigations would help in clarifying the relation of these three correlation 
measures. Along these lines, in the present paper, we have investigated the dynamics 
of quantum correlations such as quantum discord, entanglement and Bell nonlocality 
between two qubits in a lossy cavity. Very recently, entanglement, Bell nonlocality and 
quantum discord have been studied in cavity QED [HI HH1 HH EOl [221 1231 12H1 E21 [Ml 
[35l [40] . Most of them assume identical atoms, or the atom-cavity resonant case, or only 
one excitation. Here, we have investigated the effect of atom-cavity detuning and cavity 
decay rate on the corresponding correlation measures for a system which is initially 
prepared in separable or maximally entangled state. The most important results of the 
present study are: depending on the initial preparation of the system and the type of the 
detunings between the frequency of the cavity mode and the transition frequency of the 
qubits, the induced quantum discord is found to increase for a long time under cavity 
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decay loss mechanism for certain conditions, while there is no increase in entanglement 
and Bell nonlocality. In some situations, the quantum discord is induced even while the 
atoms remain separable at all times. Also, the preservation of induced entanglement 
and quantum discord was observed under strong dissipation for certain type of initial 
states. 

The paper is organized as follows. In Sec. El the model is introduced and the exact 
differential equations governing the dynamics of two atoms and single cavity mode are 
derived. The correlation measures; concurrence, CHSH-Bell inequality and quantum 
discord are also introduced, briefly. In Sees. 13.11 and I3.2[ the effect of cavity damping 
on the corresponding correlation measures are revealed for different detunings for the 
atoms initially prepared in separable or maximally entangled state. The summary of 
the important results is given in Sec. HJ 

2. The Model and Correlation Measures 

Here we consider two two-level atoms with excited states \e A ) , \e B ) and ground states 
\9a) i \9b) interacting with a single-mode cavity field of frequency u. The Hamiltonian 
of this system is given as [JT] : 

H = ^y-(\e A ) (e A \ - \g A ) (9a\) + -|^(|e B ) (e B \ - \9b) (9b\) + 

+ %gi{a \e A ) (g A \ + a ] \gA) (e A \) + hg 2 (a \e B ) (#b| + a 1 \g B ) (e B \), (1) 

where u A and u B are the transition frequencies for the atoms A and B, a and at are the 
annihilation and creation operators for the field and g\ and g 2 are the coupling constants 
to the cavity mode for the atoms A and B, respectively. 

When maximum two-photon is considered to be present in the cavity (i.e., the total 
excitation is equal to two) and if there is no dissipation in the cavity modes, the closed 
form of the total state vector at any time may be written as [41, 42J: 

|#(t)) = Ci \e A , e B , 0) + C 2 \e A , g B , 1) + C 3 \g A , e B , 1) + C 4 \g A , g B , 2) , (2) 

where Ci (i = 1,2,3,4) are the time-dependent probability amplitudes. For ex- 
ample, \g A , g B , 2) represents that the atoms are in their ground states while there 
are two photons inside the cavity. For this case, the dynamics of the atoms 
and single mode cavity field can be solved by using time-dependent Schrodinger 
equation and recently many studies have been reported for this case [HI [42} 
S3]. On the other hand, when cavity decay is taken into account, since 
the cavity decay only changes the number of photons inside the cavity, the 
atoms and the cavity mode can be found in any of the states [101 HH 14*4"] : 
\e A , e B , 0) , \e A , g B , 1) , \e A , g B , 0) , \g A , e B , 1) , \g A , e B , 0) , \g A , g B , 2) , \g A , g B , 1) , \g A , g B , 0). 
Then the density matrix of the system of two atoms interacting with a single mode cavity 
field will evolve according to the master equation [IDJ HTJ 05] : 

p = — — [H, p] io)ap — 2apo) + po)a \ , (3) 
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where a\ = |ej) (gi\ and a % _ = \gi) (e^ | (2 = A, B) are the spin inversion operators. 
Here the first term in Eq. ([3]) is generated by the Hamiltonian given in Eq. ([1]) and 
represents the coherent evolution, while the other term denotes the incoherent dynamics; 
k is the rate of loss of photons from the cavity mode. By setting {|1) = |e A ,eB,0), 
|2) = \e A ,g B ,l), |3) = \e A ,g B ,0), |4) = \g A ,e B ,l), |5) = \g A ,e B ,0), |6) = \g A ,g B ,2), 
1 7) = \g A ,g B ,l), |8) = \9a,9b,Q)} as a basis and considering the Eqs. ([1]) and ([3]), the 
differential equations governing the dynamics of two atoms and the single mode cavity 
field can be obtained easily. In our analysis we will restrict ourselves for the initial states: 
p(0) = |e A ,e fl ,0)(e A ,e B ,0|,p(0) = \e A , g B , 1) (e A , g B , 1| and p(0) = |*(0)) (*(0)|, 
where |^(0)) = \{\g A , e B ,l) + \e A , g B , 1)). One should note for these initial states 
the diagonal elements and some of the off-diagonal elements form a closed evolution. 
They can be listed as: 

Pl,l = Wl(Pl,4 - P4,l) + ^2(Pl,2 - P2,l), 

P2,2 = «V / 25 , l(p 2 ,6 - P6,2) + i>9l{.Pl,\ ~ P\,%) ~ K P2,2, 

P3,3 = igi(P3,7 ~ P7,s) + K P2,2 

p 4 ,4 = %(P4,1 - Pl,4) + «V / 25'2(P4,6 - P6,4) ~ «P4,4, 
P5,5 = W2{P5J ~ PT,5) + KP4A, 

p 6)6 = iV2gi(p 6 ,2 - P2,e) + iV2g 2 {peA - P4,e) - 2Kp 6j6 , 

P7,7 = igi(f>7fl - P3j) + m{P7,B ~ P5j) + «(2p 6 ,6 ~ P7j), 
p8,8 = ^P7,7, 

Pl,2 = - iS 2 Pl,2 + igi(V2pi,6 ~ P4,2) + m{pi,l ~ P2,2) ~ ^Pl,2, 

Pia = - ^iPi,4 + igi(pi,i - Pi,i) + ^(v^pi.e - P2a) - t^Pm, 

Pi,e = - + S 2 )pi,6 + igi(V2pi t 2 - P4,e) + ^(V^pi^ - P2,e) - «Pi,6, 
P2,4 = - - 8 2 )p2A + Wi(P2,i - V2p 6A ) + ig 2 {V2p2,6 - Pia) ~ K P2,4, 

P2,6 = - iSipzfi + iV2gi{p2,2 - P6,e) + ig2{V2p 2 A - Pi,e) - ^p 2 ,6, 

P3,5 = - i(Sl - S 2 )P3,5 - igiP7,5 + W2P3J + «/>2,4, 

P3,7 = - ^lP3,7 + «0l(P3,3 - P7,7) + W3,5 + ^ ( 2 V / 2p2,6 - jOSj), 

3 

P4,6 = - ^2P4,6 + ^l(v / 2p 4 ,2 ~ Pl,e) + ^^#2^4,4 ~ P6,e) ~ ^P4,6, 

p5,7 = - ^ 2 p 5 ,7 + igipb,3 + m{pb,5 ~ P7,j) + -{2\/2p4,6 ~ P5j), (4) 

where 5i = u A — u and <5 2 = u B — u are called the detunings and describes the frequency 
difference between the energy levels of atoms and the cavity mode, pi.j = (i \ p \ j) and 
Pi,j — Pji- We will consider two particular detuning types: for 8\ — 8\ — 8 the atoms 
will be called identical since the transition frequencies of the atoms are the same and 
the atoms will be called unidentical for 5\ = —5i = 5 case for which the individual atom 
transition frequencies are different [12]. In practice, the detunings might be slightly 
different due to broadening, in this work we will consider relatively large values of 8. 
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One should note that it is quite difficult to obtain an analytic solution for this set of 
coupled differential equations, therefore we will solve these equations numerically for 
the considered initial states. 

The density matrix of the two atoms is obtained by taking a partial trace over the 
cavity degrees of freedom: pab if) = Tr c {p} = X^=o ( n \ P \ n ) anc ^ used to calculate the 
quantum correlations as measured by CHSH-Bell inequality, concurrence and quantum 
discord between two atoms interacting with a single mode cavity field. For the considered 
initial states whose dynamics is given by Eq. (HI), the reduced density matrix of the atoms 
in the two-qubit standart basis {|1) = le^, e.e),|2) = je^, gg) , |3) = \gA, &b) , |4) = 
\qai9b)} can be calculated as: 

(a \ 
b e 



pAB(t) 



o e* c 
V d J 



(5) 



where a = p Xjl , b = p 2 ,2 + P3,3, c = p 4 ,4 + P5,5, d = p 6 ,6 + P7j + Ps,8 and e = p 2 ,4 + P3,5- 
It should be noted that the density matrix in Eq. (JSJ) has an X-form and the dynamics 
preserve its form. The considered correlation measures for this type of states can be 
calculated easily [31 [TTJ [231 12H E2] • The concurrence [26] as an entanglement measure 
is given by [XT] : 

C{t) = 2 max{0, |e| - Vad}. (6) 
The Bell nonlocality measured by the CHSH-inequality [17] can be found as [32] : 
B(t) = ma X {B 1 {t),B 2 (t)}, 

B x (t) = 2Vul + u2, B 2 (t) = 2y/ul + u3, (7) 

where ul = u3 = 4|e| 2 and u2 = (a + d — b — c) 2 . 
QD is defined as [2]: 

D(t) = I(t)-J(t), (8) 

where I(t) = S(pa) + S(ps) — S(pab) measures the total correlation between the atoms; 
S(p) = — Tr(plog 2 p) is the von-Neumann entropy and pa (pb) is the reduced density 
matrix obtained by tracing pab over the subsystem B (A). The other quantity J(t) is the 
amount of classical correlations between the atoms defined as the maximum information 
one can get about the atom A by performing measurements on the atom B, or vice versa. 
Its calculation requires maximization over one of the subsystem and generally obtaining 
an analytic result is not an easy task [3]. Recently, for the density matrix having X-form 
the explicit expressions of QD and classical correlation are reported [3j [23J, [21] . We have 
used the expressions given in Ref. [24]. The classical correlation J(t) for the density 
matrix Eq. (jSJ) is given as: 

J{t)=max{C 1 ,C 2 }, (9) 

where Cj = M(a + b) — Pj. And QD is given as: 

D(t) = win{Q 1 ,Q 2 }, (10) 



Quantum Correlations 



7 



where Qj = M(a + c) + Yh=i Kib ^°&2 KiB + Pj-> w hh \ ] AB being the eigenvalues 
of pab, Pi = M(t), P 2 = -(a log 2 a + b\og 2 b + clog 2 c + d\og 2 d) - M(a + c), 
t = (1 + ^(1 -2(c + d)) 2 + 4|e| 2 )/2 and M(a) = -alog 2 a - (1 - a) log 2 (l - a). 
We have also calculated QD numerically by considering projective measurements [2l[T4"] 
and observed that the results obtained by using Eq. (TTUT) and the numerical results agree 
perfectly. 

As is well known, concurrence is equal to 1 (0) for a maximally entangled (separable) 
state. Although QD is equal to 1 for a maximally entangled state, it may or may not be 
equal to zero for a separable state because it was shown that even some separable states 
can carry non-zero quantum discord [21 H] ■ For pure states QD is found to be equal to 
the entanglement of formation |3|. On the other hand when B(t) > 2 the CHSH-Bell 
inequality is violated which signifies that the correlations cannot be accessible by any 
classical local model [36] . 

3. Results 

3.1. Effects of Cavity Damping- Separable States 

In this section, we will examine the effect of cavity decay rate, k, on the quantum 
correlations as measured by concurrence, CHSH-Bell inequality and quantum discord. 
To do this we will solve the differential equations in Eq. fl3j) numerically for the 
initial states: p(0) = \e A ,e B ,0) (e A ,e B ,0\ , p(0) = \e A , 9b, 1} (e A , g B , 1| and p(0) = 
|tf(0)) (tf(0)|, where |tf(0)> = -±(\g A ,e B ,l) + \e A ,g B ,l)) and we will set g x = g 2 = 
g (symmetric coupling), 5\ = 5 2 = 5g for identical atoms and Si = —5 2 = 5g for 
unidentical atoms. One should note that the detunings are chosen somewhat large 
because at higher detunings the qubit-cavity energy exchange is low [UJ S2J US] and a 
high-degree of cavity- induced quantum correlations can be created [48J. 

We first consider the dynamics of concurrence and quantum discord for initially 
separable state p(0) = \e A ,e B ,0) {e A ,e B ,0\ and display C and D versus dimensionless 
time, gt in Figs. 1 and 2 for unidentical and identical atoms, respectively. For this 
initial state neither non-zero concurrence nor violation in Bell inequality ever occurs 
for unidentical qubits, so their dynamics are not plotted in Fig. 1. However, this does 
not indicate the loss of quantum correlations; QD induced by the atom-field interaction 
oscillates nearly between and 0.18 for k = (Fig. 1(a)). A similar result is found 
by Sun et al, in Ref. [12] where the authors studied the entanglement and quantum 
discord between two qubits which are independently coupled to a common environment 
modelled as an Ising spin chain. It can be seen from the analysis of Figs. 1(b)- 1(e) 
that QD is created independently of the cavity decay rate, but its maximum value 
depends on the decay rate k inversely and the effect of cavity decay on the dynamics 
of QD is to create a damped oscillatory behaviour; damping in QD is proportional to 
the cavity decay rate. Also, in order to understand the nature of non-zero quantum 
discord, we have plotted the total correlation, I(t), and classical correlation, J(t), for 
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Figure 1. Quantum discord versus gt for unidentical atoms (S± — —52 = 5g) and 
p(0) = \eA, es,0) (eA,eB,0\ initial state for the decay rates: (a) k = 0, (b) n = 0.02g, 
(c) k — 0.2g, (d) k — 2.0g and (e) k = 2Qg. Note that under these conditions 
entanglement and Bell nonlocality are not induced, thus these are not plotted here 
and the inset in (c) represents the total I(t) and classical J(t) correlations versus gt 
for k = 0.2g. 



n = 0.2g in the inset of Fig. 1(c). It should be noted that for this initial state, the total 
as well as classical correlations are also induced and I(t) is always greater than J(t) 
as a consequence non-zero quantum discord is present. This inset also demonstrates 
that I(t) and J(t) are also damped oscillatorily as D(t). Correlations in the system 
are gradually lost to the reservoir because of the cavity decay as found in the system 
studied in Ref. |18|. 
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Figure 2. Quantum discord and concurrence (insets in (a), (b) and (c)) versus gt for 
identical atoms (Si = 62 = 5g) and p(0) = |e^,eB,0) (e^,eB,0| initial state for the 
decay rates: (a) k = 0, (b) k = 0.02g, (c) k = 0.2g, (d) re = 2.0g and (e) « = 20g. 
Note that under these conditions Bell nonlocality is not induced, and for k = 2g and 
k = 20g, C(t) has no dynamics as well, thus they are not plotted here. The insets are 
plotted under the same conditions as D(t). 



We have calculated the dynamics of CHSH inequality, concurrence and quantum 
discord for the same initial state p(0) = |e^, e#, 0) (e^, e#, 0| for the identical atoms (Si = 
$2 = %) case to understand the effect of the type of detuning for the same leaky cavity 
and plot the resulting D and C as a function of gt in Figs. 2(a)-2(e) and insets of 
Figs. 2(a)-2(c), respectively. As can be seen from the insets of Figs. 2(a)-2(c), non-zero 
entanglement, although small, is created which decays to zero as cavity decay rate is 
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increased and the induced entanglement has no dynamics for higher values of k (k — 2g 
and k = 20g). Bell inequality is found to be not violated for any value of k, which is 
expected because entanglement is low and generally Bell inequality violation is shown 
to be a hallmark of high C(t) [HJ [MJ [35], so B(t) is not displayed in Fig. 2. One 
should note the difference in the dynamics of QD in the leaky cavity for unidentical and 
identical atom cases (Figs. 1 and 2). While QD has an oscillatory decay type dynamics 
for the unidentical atoms (Fig. 1), it increases until the time t = t max and then decays 
for the identical atoms. Also, for the same decay parameters, QD is found to have a 
much longer lifetime for the identical atoms. The higher leakage in the cavity decreases 
t max and the lifetime of QD after the time t = t max . On the other hand, k has no 
appreciable effect on the maximum value of D(t). Also note that the fast oscillations in 
QD are gradually smoothed for larger values of k and longer times as can be seen from 
Figs. 2(b)-2(e). 

Now we consider another initial separable state of the form p(0) = \e A , g B ) {g-Ai 9b\® 
|1) (1| which has one excitation in the field and the other excitation in one of the atoms 
to see the effect of initial partition of the excitation on the dynamics of the quantum 
correlations. The time dependent of CHSH-inequality, concurrence and quantum 
discord are displayed in the first, second and third columns of Fig. 3, respectively. 
In these figures, insets show the dynamics for unidentical atoms {5\ = —62 = 5g) 
while main subfigures are for the identical atoms (5i = 62 = 5g). Inspection of 
subfigures demonstrates the pronounced effect of qubits being identical or unidentical. 
For identical atoms, the atom- field interaction induces high values of C(t) and D(t) 
and Bell inequality is also violated at the high values of entanglement for a high- 
quality cavity (Fig. 3(a), 3(b) and 3(c)). In the case of unidentical atoms, C(t) and 
D(t) are also created, but their maximum are small; 0.018 and 0.13, respectively, 
and B(t) is not violated. When unidentical atoms interact with a low-quality cavity, 
entanglement can cease to exist very quickly and QD shows a fast oscillatory decaying 
behaviour (see the insets in Figs. 3(e)-3(r)) which is very similar to what we have 
found for p(0) = \e A ,e B ,0) (e A ,e B , 0| initial state. On the other hand, this situation is 
completely and significantly different in the case of identical atoms. Although CHSH- 
Bell inequality goes below 2 very quickly for < k < 2.0g and no longer is violated 
for k > 2.0g, the dynamics of entanglement and quantum discord are very similar 
that both approach to constant values in an oscillatory manner; the oscillations die 
down as the cavity decay rate is increased. In effect, cavity decay entangles initially 
separable state and the resulting entanglement as well as quantum discord are robust. 
In many studies, the dissipation of cavity fields were found to play a constructive role 
in the generation of quantum correlations as measured by concurrence and quantum 
discord [221 123 I2Z1 123 ESI ESI EDI- This situation can be understood by looking at 
the steady-state of the dynamics given in Eq. (JHJ). For g\ = g 2 = 1, k > and 
6i — 8 2 — 8, the p = equation has a solution independent of k and 5 and the atom- 
atom reduced density matrix is given as: p ss = \(\e A ,g B ) (e A ,g B \ + \gA,e B ) (g A ,e B \) - 
\{\zai9b) (9A,e B \ + \g A ,e B ) (e A ,g B \) + \\g A ,g B ) (g A ,g B \. The concurrence for this state 
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Figure 3. CHSH inequality (left plots), concurrence (middle plots) and quantum 
discord (right plots) versus gt for p(0) = \eA,gB,l) (e^, <7b, 1| initial state and n = 
((a), (b), (c)), K = OMg ((d), (e), (f)), K = 0.2g ((g), (h), (k)), k — 1g ((1), (m), (n)) 
and k = 20g ((o), (p), (r)). Here the main subfigures are for identical atoms with 
Si = 62 = 5g, while the insets are for unidentical atoms with Si = —S2 = 5g under the 
same decay rate parameters. Note that for unidentical atoms B(t) is not violated, so 
they are not reported here. 
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is 0.50 and quantum discord of it can be calculated from Eq. (fIU|) to be 0.412 which are 
seen clearly from Figs. 3(p) and 3(r). Also, we want to emphasize that the dynamics 
of B, C and D between identical atoms for the initial state |^/(0)) = \eA,gB, 1) is the 
same as that for |\I/(0)) = \gA, £b, 1), since we assume symmetric coupling (i.e., g\ = g 2 ) 
between the atoms and the single mode cavity field. 
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Figure 4. The populations pi,i(t) (i — 1,2,..., 8) and the absolute value of the 
coherences P2,4,(t), p3,s(i) given by Eq. (j4|) versus gt for the initial state p(0) — 
|eyi,<7B,l) (eA,9B,l\ and K = 20g for identical (S± = 62 = 5g) (a) and unidentical 
(Si = —5 2 = 5g) (b) atoms, (c) shows the atomic purity (Trp AB ) for the density 
matrix Eq. ([5]) versus gt for identical and unidentical atoms. Here the insets are for 
small time regions. 



To further elucidate the asymptotic constant entanglement and quantum discord 
for the initial state p(0) = \eA,gB, 1) (za,Qb, 1 1 , we display the time dependence of the 
density matrix elements of the atom-field system given by Eq. (pE} and the purity of the 
atomic subsystem for the density matrix Eq. (E]) for identical and unidentical detuning 
cases in Figs. 4(a)-4(c) for k = 20g. As can be seen from Fig. 4(a), all populations, except 
P3,3,P5,5 and ps,8, become zero for a short time while ps,8 goes to 1/2 and p 3;3 = p 5j5 
go to 1/4 asymptotically. The key mechanism for the creation of the entanglement and 
QD in this case seems to be the coherence generated by the interaction of the atom 
with the dissipative cavity. The absolute value of p3 t 5 also reaches to 1/4 asymptotically 
as can be seen from the Fig. 4(a). The evolved state is a trapping state which has the 
property that the probabilities of the different decay channels interfere destructively 
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and the state is robust against the cavity loss mechanism [5TI 152] . Generation or 
protecting entanglement by population trapping has been discussed by a number of 
groups before (HI HHl ETJ [52j [53] . The magnitude of k has no effect on the long time 
limit of P3,3,P5,5 and p 8j8 , it only affects how long other population components can be 
nonzero. The dynamics of atomic purity (which is equal to 1 for a pure state and 1/4 for 
a maximally mixed state in four dimensional Hilbert space) displayed in Fig. 4(c) as the 
solid line for 8\ = 62 = 5g indicates that the state approaches to a partially mixed state 
which has high entanglement and quantum discord as indicated before in Figs. 3(p) 
and 3(r). For unidentical atoms, the short time dynamics is similar to the identical 
detuning case that all states are populated for a short time, but the asymptotic state is 
now p(t —7- 00) = \g A , g B , 0) (g A , gB, 0| as can be seen from Fig. 4(b). The atomic purity 
(dotted line in Fig. 4(c)) indicates that the state of the atoms firstly goes to a mixed 
state in a short time and then approaches to the pure state p = \gA,gs) (9A,gB\ which is 
the reason of the loss of correlations between unidentical atoms in the asymptotic limit. 

3.2. Effects of Cavity Decay- Maximally Entangled Initial State 

In this subsection, we will investigate the effect of initial correlations in the atom-atom 
part of the state on the dynamics of quantum correlations. 

We take |\I/(0)) = ^(Ig^es) + \eA,gs)) ® |1) as the initial state which is the 
maximally entangled Bell state for the atomic system. QD and Bell inequality violation 
are also maximal for this initial state. We display B, C and D as a function of gt for 
several cavity decay rates for unidentical (Fig. 5) and identical (Fig. 6) atoms. At k = 0, 
all correlations show oscillatory behaviour; for the unidentical atoms B(t), C{t) and D{t) 
oscillate nearly in the ranges [2^2, 2.2], [1.0, 0.63] and [1.0, 0.58], respectively (Fig. 5). 
On the other hand, £>(£) and C(t) for k = have periodic sudden deaths and births for 
the identical detunings of the atoms (Figs. 6(a) and 6(b)), while D{t) oscillates nearly 
between 0.3 and 1 (Fig. 6(c)). For small values of k (k = 0.2g), the effect of cavity decay 
on the dynamics of correlations can be described as underdamped oscillations for both 
identical and unidentical detuning cases. As k is increased to 2g, the oscillation becomes 
overdamped and finally at large k's (k, = 20g) damping becomes critical, i.e. C(t) and 
D{f) show exponential decay for the unidentical atom cases. For the identical atoms 
at k = 2g and k = 20g concurrence has sudden death, while QD decays exponentially. 
One should note that Bell nonlocality is the most fragile quantum correlation considered 
here; it suffers sudden death in a short time as long as k > 0. 

4. Conclusion 

We have investigated the dynamics of quantum correlations as quantified by 
entanglement, Bell nonlocality and quantum discord for two qubits embedded in a 
common leaky cavity for three different initial states with double excitations. The qubits 
were considered as identical or unidentical based on the detuning of their transition 
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Figure 5. CHSH inequality (a), concurrence (b) and quantum discord (c) versus 
gt for unidentical atoms, Si = -S 2 = 5.g and p(0) = |*(0)) (*(0)| where |*(0)) = 
~^{\9a, e_B, 1) + \ca,9bA))' Here all subfigures include plots for k = 0, K = 0.2g, 
K = 2g and k — 20g. 

frequency with respect to the frequency of the field mode. The main findings of the 
study can be summarized as: Qubits being identical or not have profound influence on 
the dynamics of the considered correlations; starting from the same initial state, the 
correlations show very different behaviour depending on whether the detunings of the 
individual qubits satisfy Si = 62 = 6 or 5\ = — 62 = 5. For unidentical qubits initially in 
their excited states and interacting with a high or low quality cavity, quantum discord 
is found to be induced for all considered cavity decay rates and it decays oscillatorily 
for a leaky cavity although the qubits remain unentangled for all times. For the same 
initial state and the leaky cavity, QD is found to increase in the long time for identical 
qubits, while entanglement is either "never created" or decays to zero in a short time. 

One of the most interesting finding is the constant high entanglement and quantum 
discord state which is obtained as the long-time limit of the dynamics of p(0) = 
\ e A,9B) {^a,9b\ ® |1) (1| initial state for the identical qubit case. Here the final state is 
independent of the cavity decay rate (if k > 0) and the magnitude of the detuning and 
has high concurrence (C = 0.5) and quantum discord (D = 0.412). From an analysis 
of the time dependence of the density matrix elements it is seen that the considered 
initial state leads to a trapping state which happen to have high degree of quantum 
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correlations. On the other hand, for the same initial state, the correlations between the 
unidentical atoms are found to be destroyed quickly even for a small cavity decay rate. 

For the atomic Bell-type initial state we have found ESD for the identical qubits 
interacting with a low quality cavity, while concurrence decays exponentially for the 
unidentical qubits. Quantum discord is found to have asymptotical decay behaviour 
independent of the type of atoms for this initial state and leaky cavity. On the other 
hand, CHSH-Bell inequality violations survive only for a finite time for k > for both 
identical and unidentical qubits. It is interesting to note that steady-state nonzero 
correlation state of |^(0)) = le^Ps) <8> |1) (or |\l/(0)) = |^,e_B,l)) is not observed for 
the 1^(0)) = ^(15^4, eg) + |eA, <7b)) ® |1). This shows that the stationary states depend 
strongly on the initial conditions [49J. 

We have also demonstrated that Bell nonlocality as manifested by CHSH inequality 
is the most fragile correlation measure under dissipative environment and is violated only 
at the high values of entanglement. 
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